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We study the properties of massless and massive quarks coupled with a scalar and pseudoscalar 
boson at finite temperature in Yukawa models at the one-loop order. The behavior of the spectral 
function and the pole structure of the propagator are analyzed as functions of temperature T and 
the quark mass m/. It is shown that the three-peak structure of the spectral function found in a 
previous work for massless quarks is formed at temperatures comparable to the boson mass even 
for finite m/, but gradually ceases to exist as m/ becomes larger. We identify the three poles of 
the quark propagator corresponding to the collective excitations of the quark in the complex energy 
plane. It is shown that the three trajectories made by the poles along with a variation of T undergo 
a structural rearrangement at a critical quark mass when m/ is increased. This suggests that the 
physics content of the collective quark excitations is changed in a drastic way at this point. The 
results are nicely accounted for with the notion of the level mixing induced by a resonant scattering 
of the massive boson with quarks and holes of thermally excited anti-quarks. 
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I. INTRODUCTION 

Hadronic matter undergoes deconfinement and chi- 
ral phase transitions to the quark-gluon plasma (QGP) 
phase at finite temperature (T) . Owing to the asymptotic 
freedom of QCD, the QGP phase at asymptotically high 
T is composed of approximately free quarks and gluons, 

I where a perturbation expansion is valid. It is known that 
the leading order of such an expansion is obtained with 
the hard thermal loop (HTL) approximation [T] which 
enables us to calculate the quark and gluon propagators 
in a gauge-invariant way. One of notable features ob- 
tained in the HTL approximation is that the quarks and 
gluons have collective excitations with mass gaps called 

I thermal masses proportional to gT, where g is the gauge 
coupling [3,[1|. 

Recently, properties of the QGP phase near the crit- 
ical temperature (Tc) acquire much interest. Measure- 
ments of the elliptic flow of hadrons in the heavy-ion col- 
lisions were made at the Relativistic Heavy Ion Collider 
(RHIC), and the subsequent phenomenological analyses 
have shown that the experimental results are in remark- 
ably good agreement with the predictions of the ideal 
fluid dynamics [I]- This result suggests that the created 
matter is a strongly coupled system. The present pa- 
per is concerned with the properties of quarks in such 
a system, which constitute the basic degrees of freedom 
of the created matter in the deconfined phase together 
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with gluons. We explore possible quasi-particle picture 
of quarks in the QGP phase near Tc- 

It is not a priori clear whether quarks manifest them- 
selves as well-defined quasi-particles in such a nonper- 
turbative region; quasi-particles in Landau's sense corre- 
spond to peaks with a small width in the spectral func- 
tion with relevant quantum numbers. In the phenomeno- 
logical side, the success of the recombination model to 
describe the RHIC data Q tells us that the quark quasi- 
particles play an important role in some stage of the 
evolution of the created matter in RHIC. In the the- 
oretical side, a recent lattice QCD simulation 6] indi- 
cates the existence of the quasi-particles of quarks with 
small decay width, though in quenched approximation. 
Quasi-particle properties of quarks are also obtained in 
a strong coupling gauge theory based on the Schwinger- 
Dyson equation at finite T, unless the coupling is too 
large The quarks thus will exist as one of the basic 
degrees of freedom of the system even near Tc, and hence 
revealing their nature will provide us with important in- 
formation necessary to understand the nature of QGP 
phase at the strong coupling. 

To study the quasi-particle properties of quarks around 
Tc, we notice the possible existence of hadronic excita- 
tions even in the QGP phase near Tc [1, [§]. The nu- 
merical results in lattice QCD suggest the existence of 
such states in the heavy-quark sector ■ The existence 
of the mesonic excitations in the light-quark sector was 
also predicted as being the soft modes associated with 
the chiral transition [8]. If such bosonic modes are light 
enough so that many particles are thermally excited in 
the system, they can in turn strongly affect excitation 
spectra of quarks. Part of the authors in the present pa- 
per (M.K., T.K. and Y.N.) considered such possibility 
in Refs. [ll|, ■ They first evaluated the effect of the 
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chiral soft modes on quarks in a chiral effective model 
[m where the soft modes are described as composite 
bosons composed of a quark and anti-quark: An interest- 
ing result obtained there was that there appears a novel 
three-peak structure in the quark spectral function near 
but above Tc. It was also pointed out that the coupling 
of quarks with bosonic excitations which have a finite 
mass and small width is essential for the formation of 
the multi-peak structure. In fact, it was later shown 
that the three-peak structure in the quark spectral func- 
tion emerges at intermediate temperatures even in the 
Yukawa models composed of a massless fermion and a 
massive boson with a zero width, irrespective whether 
the boson is of scalar, pseudoscalar, vector or axial- vector 
type. 

In the previous work [ll|, the massless quarks 
were exclusively adopted to explore the effects of the soft 
modes on the quarks in the most ideal case where the 
chiral symmetry is exact. In reality, the quarks in the 
QGP phase near Tc should have current or dynamical 
Dirac masses which may be as large as the order of 
and may not be negligible. In this paper, we examine 
how the finite Dirac mass of the quark modifies the spec- 
tral properties, employing the Yukawa model composed 
of a massive quark with a small mass and a massive bo- 
son. We shall show that the three-peak structure found 
in Ref. 12] survives even with finite mf, although one of 
the three peaks is gradually suppressed as mf is larger. 

In addition to the analysis of the spectral function, we 
newly explore the pole structure of the quark propaga- 
tor in the complex-energy plane and their residues. We 
show that there always exist three poles corresponding to 
the peaks in the spectral function in the complex-energy 
plane at finite T. It is found that the trajectories of these 
poles along with an increase of T undergoes a structural 
change at a critical value of m/ as mf is increased. We 
give a physical interpretation for these properties of the 
quark spectrum at finite mf adapting the notion of the 
level repulsion in quantum mechanics. 

This paper is organized as follows. In Sees. [TIllIVl we 
study the spectral properties of the quark coupled with 
a scalar boson. After introducing the Yukawa model in 
Sec. mi we present the numerical result for the quark 
spectral function in Sec. IIIII and the pole structure of 
the propagator in Sec. IIVI We then discuss the spec- 
tral properties of the quark coupled with a pseudoscalar 
boson in Sec. El In Sec. IVI[ we elucidate the physical 
mechanism of the m /-dependence of the quark spectrum. 
Section [VIII is devoted to the summary and some discus- 
sions. In Appendix [XI we present the renornialization 
procedure of the quark propagator using the subtracted 
dispersion relations. In Appendix [Bl we introduce an an- 
alytic toy model and elucidate that a two-peak structure 
of the quark self-energy is essential for the appearance of 
the three-peak structure in the quark spectrum. To sim- 
plify the analysis, we restrict the following calculations 
to zero momentum to concentrate on the T dependence. 
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FIG. 1: The one-loop quark self energy. The solid and dashed 
lines represent the quark and the scalar boson, respectively. 4- 
momenta P and K denote P = (p = 0, oj^) and K — {k,u>n), 
respectively. 



II. YUKAWA MODEL COUPLED WITH A 
SCALAR BOSON 

In this section, we formulate a Yukawa model with a 
massive quark and a massive scalar boson and evaluate 
the quark self-energy at one loop. The quark spectral 
function and some other quantities needed in the follow- 
ing sections are also introduced. 

We start from the following Lagrangian composed of a 
quark and a scalar boson, 

C = - mf)^ + i - ml] - g^^, (2.1) 

with the quark mass to/, the boson mass to;, and the 
coupling constant g. If the boson is a pseudoscalar, g 
should be replaced with i^^g. This replacement may lead 
to non-trivial differences because the quarks are massive. 
The pseudoscalar case will be analyzed in Sec. |Vl 

The spectral properties of the quark at zero momentum 
can be extracted from the retarded quark propagator 

G^{io) = [{io + Z7?) - TO/ - S« (c^)] , (2.2) 

with the retarded self-energy I]^ (w). We calculate the 
quark self-energy at finite T at one loop, which is dia- 
grammatically shown in Fig. [1] In the imaginary-time 
formalism, it is given by 

f.{iuj„i) = -g^T^^ I go(k, ia;„)2?(-k, i^m-i^n), 

(2.3) 

where 00 (k, itt^ri) — [i'j-'„7° — k-7 — to/] ^ andX'(k, w„) = 
[{iVnY — — the Matsubara Green functions 

for the free quark and the free scalar boson, respectively, 
and ujn — (2rt + 1)itT and Un — 2m:T are the Matsubara 
frequencies for fermions and bosons, respectively. Taking 
the summation over the Matsubara frequency and mak- 
ing the analytic continuation iujn uj + irj, one obtains 
the retarded self-energy I]^(w) = ^{ii^n)\iuj„=Lu+iri- 

The self-energy E^(a;) has an ultraviolet divergence 
which originates from the T- independent part Ej:^q(w) = 
\imrp^Q+ Y^^{uj). We remove this divergence by impos- 
ing the on-shell renormalization condition on the T- 
independent part of the quark propagator; for the de- 
tails, see Appendix A. The T-dependent part, Sj;_^q((jj) = 
I]^{uj) — Ey^Q(a;), on the other hand, is free from diver- 
gence [l2| . For the numerical calculation, it is convenient 
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FIG. 2: The decay processes corresponding to the terms (II) 
and (III) in Eq. ((IS}. 
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FIG. 3: The supports of the terms (I)- (IV) in Eq. ((gS}. Here, 
M+ = mb + mf and Af_ = \mb ~ rnf\ are boundaries of the 
supports. 

to first calculate ImEy_^Q(tj) and then determine the real 
part via the Kramers-Kronig relation 



dz Inil]^_^g(z) 



(2.4) 



where V denotes the principal value. 

The quark propagator for zero momentum is decom- 
posed in the following way, 

G^iuj) = G+(w)A+7" + G-(t^)A_7°, (2.5) 

with projection operators A± = (1 ± 7°)/2 onto spinors 
whose chirality is equal(-l-) or opposite(— ) to the helicity. 
We call +(— )-sector as 'quark'('anti-quark') sector. Here 

G±{cj) = iTr [G^(a;)7°A±] 

= [cj + i?7T"i/ -I]±(w)]~\ (2.6) 

with T,±{u)) = tr[S^(a;)7°A±]. From the charge conju- 
gation symmetry, we have a relation between G± 



G+{uj) = -Gl(-tj). 



(2.7) 



Because of this relation, we consider only the 'quark' 
sector G+{uj). For nif — 0, the propagator satisfies 
G±{uj) = —G'^{—uj). Notice that at finite quark num- 
ber density, the charge conjugation symmetry is violated 
and Eq. (|2.7p is no longer valid. 

The imaginary part of S_|-(a-') is given by 

Tr[A+A+(k)](5(w - Ef ^ E,)[l + n - /] (I) 
+TT[K+K_{\<i)]5{u: + Ef-Eb)[n + J] (II) 
+Tr[A+A+(k)]<5(w - Ef + Ei,)[n + /] (III) 
+Tr[A+A^(k)]5(co + Ef + Eb)\l + n- /]} (IV), 

(2.8) 



where Ef = [m) + k2]i/2^ = [ml + ]^Y'^, A±(k) = 
(i?/ ±7°(7 • k + m/))/(2i?/)^and / = (exp(£;//r) -f 1)^1 
and n — (exp(i?{,/T) — 1) ^ are the Fermi-Dirac and 
Bose-Einstein distribution functions, respectively. Equa- 
tion (|2.8p is proportional to the difference between the 
decay and creation rates of the quasi-quark. One can 
give a diagrammatical interpretation to each term (I)- 
(IV) [H, [l^. In Fig. [2l we show the decay processes 
included in the terms (II) and (III), which are called the 
Landau damping. The term (II) corresponds to a pair 
annihilation process between the quasi-quark (Q) and a 
thermally excited anti-quark {q) with an emission of a 
boson (6), and its inverse process. The term (III) cor- 
responds to a scattering process of the quasi-quark by a 
thermally excited boson. The Landau damping vanishes 
at T = 0, as it involves thermally excited particles in the 
initial states. As we will see later, the Landau damping 
is enhanced rapidly as T goes high and plays a significant 
role for the peak structure in the quark spectrum. 

Each term (I)-(IV) has a support on the u) axis as 
shown in Fig. [3] owing to the energy-momentum conser- 
vation. While any value of uj belongs to one support for 
TO/ = 0, the ranges Ito^ — to/| < \lo\ < mh + mf are for- 
bidden kinematically for ruf ^ Q and hence ImE-|_(a;) = 
there. In particular, the supports for the Landau damp- 
ing (II) and (III) vanish for TO/ = TOf,. 

The integral in Eq. (|2.8p can be performed analytically 
to give 



ImE+(ci;) 



.g2 {u + M+){uj - M_) 
647r uj"^ 



coth 



2 -Af2)(t^2 _^^2 
2 ' M-M. 



tanh 



cj2 - A/_M, 



+ 



4ujT 



AujT 

X [6i(cj2 - Ml) - 0{m1 - ujI] , (2.9) 

where = TO/-f toj,, M_ = |TOfc — to/|, and 0{u!) is the 
step function. 

In the following sections, we investigate the spectral 
properties of the quark by analyzing the spectral func- 
tion and poles of the quark propagator. The spectral 
function is given by p{u!) = — (l/7r)ImG^ (w) which is de- 
composed, similarly to Eq. p.Sp . as p{uj) — p+{uj)A+j^ + 
p_(w)A_7°, with 



-ilmG±(w) 

TT 



IT [uj^mf ~ Rel]±(tj))2 + Imi;±(a;)2 ' 

(2.10) 

As directly derived from Eq. (|2.7p . the charge conjugation 
symmetry leads to 



p+{uj) = p-i-uj). 



(2.11) 



In general, p+{lo) has several peaks corresponding to the 
quasi-quark excitations. Such excitations can also be ex- 
pressed by complex poles of G+(aj) which are given by 
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solving 



[G+iz)]- 



S+(z) = 0, 



(2.12) 



with a complex number z. The causality requires 
that the solutions of Eq. (|2.12p are in the lower-half 
complex-energy plane. The retarded self-energy ^+(2:) 
in Eq. (|2.12p for a lower-half complex-energy is given by 



2t [lm^+{u;)l 



(2.13) 



Notice that the analytic continuation of the retarded 
function to the lower-half plane requires the second term 
in Eq. (|2.13p to compensate a cut of the first term on 
the real axis; for details of this analytic continuation, see 
Ref. [l3| • The residue of each pole reads 



Z = 



1 - 



as+(z) 



(2.14) 



at the position of the pole z. We remark that the residue 
of a complex pole is a complex number in general. 

To understand how the peak structure in p+{uj) 
emerges, it is also convenient to use the notion of "quasi- 
pole" which is defined as zero of the real part of the in- 
verse propagator; 



oj - 



ruf - Rel]+(w) = 0. 



(2.15) 



Indeed, one sees from Eq. (|2.10p that p+(w) becomes 
large at oj satisfying Eq. (|2.15p . provided that ImS+(a;) 
is small there. 



III. THE SPECTRAL FUNCTIONS OF QUARKS 
COUPLED WITH A SCALAR BOSON 

In this section, we present the numerical results of 
the quark spectral function in the Yukawa model with 
a scalar boson. In order to see the effect of the quark 
mass ruf, we first recapitulate the result for ruf — 
[l^ . focusing on its T-dependence at zero momentum. 
We then show the numerical results for massive quarks. 
The behavior of the self-energy is also analyzed for both 
cases. We fix the coupling constant g unity throughout 
this and subsequent sections. We have checked that our 
results do not change qualitatively with the variation of 
g over a rather wide range . 



A. Quark spectrum with mf — 

In the left upper panel of Fig. IH we show the T- 
dependence of the spectral function p+(w) with mf = 0. 
p+{i-u) takes on a peak with a 5- function shape at the 
origin independently of T. In this figure, the height of 
this peak denotes Z x 10 with Z being the residue of the 
corresponding pole in the quark propagator. At T = 0, 
the residue is unity owing to the on-shell renormalization 



condition. As T is raised, the magnitude of the present 
residue decreases while two bumps newly appear at finite 
energies and eventually turn to narrow peaks. At inter- 
mediate temperatures, i.e., t = T/mt, = 1-1.5, p+{ll)) 
thus has a three-peak structure [Tll.[l^. As will be shown 
in the next section, there exist poles corresponding to the 
new two peaks, and the residues of the three poles tend 
to have almost the same magnitude in this range of T, 
which indicates that the three peaks are equally signifi- 
cant. As T is raised further, the residue of the pole at the 
origin damps quickly and p+ (w) is dominated by the two 
peaks at finite energies: These two peaks can be identi- 
fied with the normal quasi-quark and the anti-plasmino 
excitations [l3| which are familiar collective fermion ex- 
citations obtained in the HTL approximation in gauge 
theories [3]. The positions of these peaks correspond to 
the thermal masses. 

To understand the T-dependence of ^+(0-"), we ana- 
lyze the imaginary and real parts of 5]+(a;), which are 
shown in Fig. [5l For nif — 0, the self-energy S+(w) 
has symmetry properties, Iml]_|_(a;) = ImS+(— tj) and 
ReE+(w) = -ReT.+ {-uj). At T = 0, Imi;+(tj) takes 
non-zero values only for \uj\ > mi, corresponding to the 
supports of the terms (I) and (IV) in Eq. dM]). At T ^ 0, 
ImE+(a;) gets to have a support in the region \uj\ < mi, 
owing to the Landau damping, i.e. the terms (II) and 
(III), and forms two peaks there. The heights of these 
peaks grow up quickly as T is raised. 

Then owing to the Kramers-Kronig relation between 
the real and imaginary parts of the self-energy Eq. (|2.4p . 
ReS+(aj) inevitably shows a steep rise in the two regions 
of Lu corresponding to the peak regions of ImE-|_(w) for 
\uj\ < mi, [12]; their steepness and the magnitudes are 
amplified as T is raised. 

This oscillatory behavior of ReS_|_(w) affects the num- 
ber of quasi-poles. In the lower panel of Fig. [SJ we show 
a line ReS+(a;) = uj to find the solutions of Eq. (|2.15p : 
crossing points between this line and ReS+(w) give the 
quasi-poles. The figure shows that at low T, there exists 
only one quasi-pole at the origin. As T is raised, the oscil- 
latory behavior of Rel]-|_(ti;) is so enhanced that four new 
quasi-poles get to exist. However, Iml]_|_(aj) has so large 
values at two of the quasi-poles that any peak structure 
in p+(w) is not formed at these points; thus we end up 
with only the three peaks in the spectral function at the 
origin and the other two quasi-poles where Im£_|_(a;) is 
small. The moral which one should learn from the above 
discussion is that the Landau damping which causes two 
peaks in ImE+(ti;) is essential for the formation of the 
three-peak structure in p+(a;); see Appendix [Bl for an 
elucidating discussion with use of an analytic toy model 
for the self-energy. 



B. Quark spectral function with small quark mass 

In this subsection, we investigate how a finite but small 
quark mass affects and modifies the quark spectral func- 
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FIG. 4: The quark spectral functions mbp+{ijj) for mf/mb = 0,0.1,0.2 and 0.3. The heights of the peaks of the 5-function 
shape represent Z x 10 with Z being the residue. 

In the lower-left panel of Fig. [4l we show the spectral 
functions for mj/mf, = 0.2. As mentioned above, there 
exists a peak of the ^-function type at u) = mj atT = 0. 
As T is raised, this peak becomes broader and eventually 
splits into two peaks; one of which approaches the origin 
while the other connects to the normal quasi-quark exci- 
tation having a non-zero thermal mass at high T . There 
also appears a broad bump in the negative-energy re- 
gion, which only gradually gets developed into a peak at 
higher T than that in the positive-energy region. Thus, 
the number of clear peaks in at a certain T is two 

at most, and a clear three-peak structure is hardly seen. 



tion p^{ljj). 

We show p+{uj) with mf/mi, ~ 0.1 in the upper-right 
panel of Fig. ID At T = 0, p+{uj) takes on a peak with 
a (5- function shape at lu = mf; this pole always exists in 
p+{u!) at T = irrespective of the value oi mf. As T is 
raised, this peak gets to have a width and moves toward 
the origin. There also appears a shoulder in the positive- 
energy region larger than m/, which turns to a narrow 
peak at T = T/mb ^ 1.5. In the negative-energy region, 
there appears also another peak whose development is 
slower than that in the positive energy region. The three- 
peak structure is then barely seen in at r ~ 1. 
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for mf/mb = 0.2. Comparing this figure with 
Fig. O one finds that the behavior of S+(aj) hardly 
changes with this shght change of the quark mass. Since 
the quasi-poles are given by solving ReS+(aj) = co — mf , 
the dominant effect of mf thus comes from the quark 
mass term in the r.h.s. of the above equation. The 
line to determine the solutions of this equation is shifted 
lower with ruf ^ 0, as shown in Fig. [51 Owing to this 
shift, there arises an asymmetry in the positions of quasi- 
poles in positive- and negative-energy regions. The two 
quasi-poles in the positive-energy region emerge at lower 
temperature, while the quasi-pole in the negative-energy 
region at higher temperature. 



FIG. 5: The quark self-energy E+(cij) with m/ = for 
T/rub ~ 0, 0.5, 1.0 and 1.5. The upper and lower panels show 
the imaginary and real parts, respectively. 
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FIG. 6: The quark self-energy E-i-(tj) with mf/mb = 0.2 for 
several values of T. 



In the lower-right panel of Fig.[4l the spectral function 
p-i- for mf/mb = 0.3 is shown. We see that the peak 
at Lu — mf at T — smoothly connects to the normal 
quasi-quark excitation in the high T limit, which is quite 
different from the cases with mf/mb £ 0.2. Further- 
more, the three-peak structure in p+{(-u) does not appear 
for any T; this is because the peak in the negative-energy 
region develops only very slowly with increasing T. As 
we will see later in Sec. IIV[ this behavior reflects in the 
T-dependence of the poles of the quark propagator. This 
qualitative change of the T-dependence of the pole be- 
haviors can be understood in terms of the level crossing 
as will be discussed in Sec. IVII 

The effects of on the behavior of P+{lo) can be un- 
derstood more precisely by analyzing the m/-dependence 
of the self-energy and quasi-poles. In Fig. [S] we show 



IV. THE POLE STRUCTURE OF QUARK 
COUPLED WITH SCALAR BOSON 

In this section, we investigate the pole structure of 
the quark propagator G+{uj). We find the three poles 
which correspond to the peaks in p+{lo) discussed in the 
previous section, and analyze their T-dependence. It is 
found that the T-dependence of the poles shows a drastic 
change at mf/mb — 0.21, which means that the physics 
contents of the collective quark excitations are changed 
at this point. 



A. Pole structure of massless quark 

We first examine the case with m/ = 0. The poles 
of the quark propagator are given as complex roots of 

Eq. ^rm . 

We first notice that G-\-{z) has an infinite number of 
poles in the lower-half complex-energy plane. In partic- 
ular, there are two series of infinite poles near the imag- 
inary axis whose density increases infinitely as they ap- 
proach the origin. The manifestation of these poles is 
mathematically attributed to the hyperbolic cotangent 
and tangent in Eq. (|2.9p . 

where we have set mf — 0. For mb ^ 0, the argu- 
ments of them diverge as ~ 1/z at the origin, and these 
terms lead to a rapid oscillation of in the complex- 

energy plane. On the imaginary axis, for example, this 
term behaves as cot(l/C) -|-tan(l/C) near the origin with 
^ = —iz, and thus ImE-|_(C) has infinite number of diver- 
gences around the origin. This behavior also brings about 
singularities of ReS+(C) through the Kramers-Kronig re- 
lation Eq. (|2.4p , and leads to the infinite number of roots 
of Eq. (|2.12p near the imaginary axis. These poles, how- 
ever, do not correspond to any peaks in the spectral func- 
tion, and thus has no physical significance. Because this 
oscillation comes from the T-dependent term, this is spe- 
cific at finite T. 
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FIG. 7: The T-dependence of the poles (A)-(C) with m/ = 0. 
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FIG. 8: The T-dependence of ReZ of the poles (A)-(C) with 
nif — 0. The summation of KeZ, Eq. (|4.4|l . is also shown. 



Besides these poles, we have found that there exist 
isolated three poles which have a physical significance in 
the sense that they all correspond to the peaks in p^{uj). 
One of them is that existing at the origin independently 
of T. We label this pole as (A) in the following. For r > 
2.0, we find the other two poles associated with the sharp 
peaks in p^{u;); one has a positive real part and is labeled 
as (B), while the other has a negative real part and is 
labeled as (C). As T is lowered, the poles (B) and (C) 
move toward the region containing the above-mentioned 
infinite poles and become hard to distinguish from them 
numerically: For r < 0.4, our numerical algorithm failed 
to identify these two poles, and the pole search is stopped 
at T = 0.4. In the following, we concentrate on the poles 
(A)-(C), since they have physical significance: Indeed it 
will be shown that the peak structure of G+{u!) is well 
described solely by these three poles. 

In Fig. [71 we show the T-dependence of the positions 



of the poles (A)-(C). The real parts of the residues at 
these poles, KeZA, RgZb and KeZc, are also shown in 
Fig. [S] One sees from these figures that for relatively 
low T, |Imzs| and jlmzd are large and KeZs and HeZc 
are small, which is consistent with the result in the pre- 
vious section that p+{lli) is dominated by the pole (A) 
at sufficiently low T. For r ~ 1.2, all the residues have 
the same value, which indicates that the collective exci- 
tations corresponding to these three poles have a similar 
strength. For high temperatures, llmz^l and jlmzd be- 
come smaller in accordance with the formation of the 
sharp peaks in p+{uj). 

In order to check that the poles (A)-(C) certainly cor- 
respond to the peaks in we compare the spec- 
tral function with that constructed from a Breit-Wigner 
(BW) approximation 



TT ^ — ' UJ — Zj 

i=A,B,C 



(4.2) 



This approximated and the original spectral functions 
for T — 0.5, 1 and 1.5 are plotted in the upper panels of 
Fig.[9l One sees that p^{uj) well reproduces the original 
spectral function p+{uj), especially for higher T. We thus 
conclude that the three poles well represent the excitation 
property of the quasi-quark. The BW approximation, 
however, fails to reproduce the precise shape of p+{uj) 
around the origin, which is due to the contributions of 
the infinite series of poles near the imaginary axis. 
The fermion spectral function satisfies the sum rule 



rfcj/9+(cij) = 1. 



(4.3) 



If the poles (A)-(C) well describe the whole analytic 
structure of the propagator, the sum of the residues 
^ ■ Zi should take a value close to unity, provided that 
a possible violation of the sum rule owing to a renormal- 
ization procedure is neglig ible 

Here, we first notice that, with the BW approximation 
Eq. ((42)) . the left hand side of Eq. (|43|) is given by, 



dujpl'^iLj) = 



E 



Rc Zi = Zto 



(4.4) 



which tells us that only the real part of the residue con- 
tributes to the strength of the spectrum. The T depen- 
dence of Zxot as well as ReZi {i =A, B and C) are shown 
in Fig. [S] One sees that the sum Ztoi gives a value close 
to unity in the whole range of T with a small deviation, 
about 20% at most. We remark that if the values of 
residues of the infinite series of poles around the imagi- 
nary axis is added to Zxot the deviation from unity be- 
comes much smaller. 

Finally, we examine the poles (A)-(C) for t ^ 1. In 
Fig. [lOl we show the poles (A)-(C) up to r ~ 40. The 
poles (B) and (C) reach the real axis at w = mi, when 
T ~ 4.2. This is due to a suppression of the phase space of 
the decay processes at w = mi,. For higher T, |Imz/rnb| 
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FIG. 9: The comparison of the spectral functions with the Breit-Wigner approximations {(^) and the original one p+{u 
for several quark masses and temperatures. 
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FIG. 10: The T-dependence of the poles (A)-(C) with mf = 
up to r ~ 40. The points denote the position of poles for 
r = 5, 10, 15, 



for the poles (B) and (C) increase monotonically with T. 
At sufficiently high T, the real parts of the poles (B) and 
(C) take a value close to those determined in the HTL 
approximation having the thermal mass ttit — gT /A p^ . 
although there is a small deviation owing to the strong 
coupling 3=1 taken here. The behavior of the imaginary 
part of each pole is also consistent with that in the HTL 
approximation which gives the width of order grriT 3- 
Although there always exists the pole (A) at the origin, 
the residue of this pole decreases and becomes negligible 
in the high-T limit. 



B. Pole structure of massive quark 

In this subsection, we examine the T-dependence of the 
pole structure of the propagator for the massive quarks 
by varying the quark mass (m^/mf, = 0.1 ^ 0.3). We 
shall show how the small but non- vanishing quark mass 
significantly affects the physics contents of the quasi- 
quark excitations at finite T . 

First of all, we have found that there exist the three 
poles of physical significance which correspond to a peak 
or bump in the spectral function even for the mas- 
sive quarks, too: The lower panels of Fig. [9l show the 
spectral functions in the BW approximation {lu) for 
mf/nii, = 0.1, 0.2 and 0.3 with fixed r = 1. We see that 
the BW approximation solely with the three poles well re- 
produces the original spectral functions even for the cases 
of massive quarks, which means that the three poles well 
represent the physical excitations of quarks even with fi- 
nite ruf. Here we mention that an infinite number of 
poles due to a kinematical origin exist also for finite to/ 
around the imaginary axis as in the case of m/ = 0. 

We show the T dependence of the poles for m y /mf, = 
0.1 in Fig.[TT] The figure shows that the pole (A) moves 
toward the origin as T is raised, while the poles (B) and 
(C) move toward the real axis almost symmetrically with 
respect to the imaginary axis: The small asymmetry of 
the two trajectories is caused by the fact that |Im0s| 
is always smaller than |Imzc| in the range of t shown 
in the figure. This inequality for the imaginary parts 
is consistent with the behavior of p+{oj) that the peak 
in the positive-energy region is sharper than that in the 
negative-energy one. In Fig. [T^ we show KeZi of each 
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FIG. 11: The T-dependence of the poles (A)~(C) with m/ 
0.1. 



FIG. 13: The T-dependence of the poles (A)-(C) with 
mf /mb = 0.2. 
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FIG. 12: The T-dependence of ReZ with ruf = 0.1. 



pole in the case of nif/nih — 0.1. The T-dependence of 
each residue is qualitatively the same as that for m/ = 0. 
We notice that the residues for all the poles have similar 
values at T ~ 1 where a three-peak structure is barely 
seen in p^(uj). For temperatures much higher than those 
shown in Fig. [Til the poles (B) and (C) eventually reach 
the real axis at \KezB,c\ — "nib — mj- As T rises further, 
they move on the real axis until jRezB.cl — "nib -I- m/ at 
which the poles start to leave the real axis. This move- 
ment of the poles along the real axis for some temper- 
atures reflects the support structure of the ImS+(cL;) as 
shown in Fig. [31 For r 1, the positions of the poles 
approach those for mf = shown in Fig. 1101 

Let us turn to the case with mf/mh — 0.2. The T- 
dependence of the poles in this case is shown in Figs. [121 
Although the topology of the set of the trajectories and 
the behaviors of the poles along with the variation of 
T are qualitatively the same as before, a notable point 
is that the trajectories of the poles (A) and (B) are so 



FIG. 14: The T-dependence of ReZ with mf/nib = 0.2. 



deformed that they tend to get close to each other around 
T = 0.6. As we will see later, this is actually a precursory 
deformation leading to a structural rearrangement of the 
trajectories to be seen for the larger mf/rrih. 

The T-dependence of the pole residues ReZi [i — A, B, 
C) is shown in Fig.[T31 We see that the behavior of them 
are quite different from those for smaller quark masses: 
As T is raised, HcZa first increases and then rapidly 
drops at r ~ 0.6, while KeZs behaves oppositely. The 
sum of the residues Zxot is, however, almost unity for any 
T. It should be noted that although ReZs is negative 
around r = 0.6, the pole (B) does not correspond to 
any peak in the spectral function nor hence describe a 
physical excitation. 

Finally let us examine the case with mf/nib = 0.3. 
Fig. [151 shows how the three poles move along with the 
increase of T. At a first glance, one can see that a dras- 
tic change has occurred by a small increase of mf/mb 
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FIG. 15: The T-dependence of the poles (A)-(C) with 
mf/rrib — 0.3. 
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FIG. 16: The T-dependence of ReZ of the poles (A)-(C) 
with mf/nib = 0.3. 



from 0.2 in the topology of the set of the three trajec- 
tories as well as in the T-dependence of the poles. The 
pole (A), which is identified as one smoothly connecting 
to cj = TO^ at T = 0, moves with increasing T toward 
large positive energies with a small or vanishing imagi- 
nary part, as pole (B) did in the previous cases, instead 
of going to the origin as the pole (A) did before. On the 
other hand, the pole (B) first appears with a large imag- 
inary part near the imaginary axis at r = 0.4 and tends 
to move toward the real axis at high temperatures as in 
the cases of irif /mt, < 0.2, but then changes the direction 
and approaches the origin as T is raised further, as the 
pole (A) did in the previous cases. We now see that the 
poles (A) and (B) at high temperatures has exchanged 
their characters with each other; this exchange should 
have occurred at a critical quark mass nifc which is lo- 
cated between 0.2mb and O.Smf,. The character exchange 



of the poles is a kind of level crossing familiar in Quantum 
mechanics. In geometrical terms, when m/ = m/c, the 
trajectories of pole (A) and (B) are so deformed that the 
two trajectories would touch each other at some temper- 
ature (r ~ 0.6), and then for larger quark masses than 
nifc, the two trajectories would be rearranged to form 
trajectories as shown in Fig. 1151 Indeed, we have numer- 
ically confirmed that the trajectories of the poles (A) and 
(B) get close to and then touch each other eventually at 
mf/mb ~ 0.21, i.e., nifc — 0.2mb. The T-dependence of 
KeZ of each pole for mj /mi, — 0.3, shown in Fig. [111 can 
also be understood naturally with the notion of the level 
crossing. 



V. SPECTRAL PROPERTIES OF QUARKS 
COUPLED WITH PSEUDOSCALAR BOSONS 

In this section, we replace the scalar boson with a 
pseudoscalar one and investigate how the behavior of 
the spectral function and the pole structure change. 
Throughout this section, we put a superscript "PS" to all 
the functions where PS is for pseudoscalar boson. When 
we refer to the functions calculated in the previous sec- 
tions, we put a superscript "S" where S is for scalar bo- 
son. 

The Lagrangian density of the Yukawa model com- 
posed of a quark and a pseudoscalar boson is given by, 
1 



(5-1) 

The quark self-energy in this model in the imaginary time 
formalism at one loop reads 



175^0 (k,iw„) 



(27r)3 

X i752?(-k, ium ~ iuJn)- (5.2) 



The difference between Eqs. (|5.2p and (|2.3p is the two 
factors of 175 only. As one can easily check, this difference 
leads to the following simple relations for the retarded 
self-energies. 



ReEf (w) 
ImI]PS(a;) 



-ReS'=_(-w), 



(5.3) 
(5.4) 



For irif — 0, one immediately obtains Sl^^(tj) — I]^(tj), 
and thereby the type of bosons does not affect the quark 
propagator at all [12| in accordance with the chiral sym- 
metry. 

In Fig. 1171 we show the imaginary and real parts of the 
self-energy with the scalar and the pseudoscalar bosons 
for mf/nib — 0.3 and r = 1.5. We see that the behavior 
of them is qualitatively the same as that for the scalar 
boson. 

From this feature, one naturally expects that the re- 
sulting spectral function and the pole structure are also 
insensitive to the type of the boson. In Fig. [THl we show 
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FIG. 17: The imaginary and real parts of the self-energies 
E'^(w) and E^^(cj) for nif/mt = 0.3 at r = 1.5. The dashed 
line is the line uj ~ nif. 
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FIG. 18: The comparison of the spectral functions of the 
quark coupled with the scalar and pseudoscalar bosons. The 
labels "S" and "PS" stand for the scalar and pseudoscalar 
bosons, respectively. 



the spectral functions p^^{uj) and jO^(aj) for mf/nih = 0.3 
and T — 1.5. One sees that the peak structures of these 
spectral functions are qualitatively the same, indeed. 
The pole position is shown in Fig. [12] for mf/rrib — 0.2 
and 0.3. We finds again that each pole behaves simi- 
larly as before. Quantitatively, |Imz^| is always smaller 
than that obtained in the previous section, while jlmzsj 
is larger. The level crossing between the poles (A) and 
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FIG. 19: The T-dependence of the poles. The upper and 
lower panels are for mj /rub = 0.2 and 0.3, respectively. 

(B) occurs for this case also at nif/mi, ~ 0.23, which is 
slightly larger than that for the scalar boson. 

VI. DISCUSSIONS; LEVEL MIXING INDUCED 
BY RESONANT SCATTERING 

In the preceding sections, we have studied the spectral 
properties of quarks coupled with scalar and pseudoscalar 
bosons by examining the spectral function p+{uj) and the 
pole structure of the quark propagator. We have found 
that, although the three-peak structure for ruf = sur- 
vives even for finite nif, such structure gradually ceases 
to exist and the number of clear poles in p+{uj) becomes 
two at most for rrif /mb > 0.2. The behavior of p+{w) and 
the position of poles at finite m / show a drastic change 
as if there is a level crossing at mf/mt ~ 0.21 (0.23) 
when coupled with a scalar (pseudoscalar) boson. In this 
section, we elucidate the mechanism for realizing this be- 
havior in terms of the notion of the so called resonant 
scattering [15]. 

We first recall that the formation of two peaks in 
lmY]^{uj) plays a decisive role to cause the three-peak 
structure in p+{uj) (See, Sec. IIIII and App. |B]). The de- 
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FIG. 20: The energy levels for mf/mt = 0.2 (left) and 0.3 
(right). The solid lines are Lj(")(ifc) for w > and t^("^'(fc) 
for uj < 0. The crosses indicate the free-quark energy at zero 
momentum, u = rrif. The dashed lines e'^^ and £^™' de- 
note the most probable energy levels of uj^^^^k) and ijj^™\k), 
respectively.. 



cay rates corresponding to these two peaks are due to the 
Landau damping, which is diagrammatically depicted in 
Fig. [21 The term (II) is the pair annihilation process of 
the quasi-quark (Q) and a thermally excited anti-quark 
(g) into a boson b, Q + q ^ b, and its inverse process. 
Here, we notice that the annihilation of a thermally ex- 
cited anti-quark can be regarded as the creation of a hole 
with a positive quark number in the anti-quark distribu- 
tion. With this interpretation, the term (II) is schemati- 
cally described as Q — * ^hoie ^ ^. ^^^^ quasi-quark is con- 
verted to a hole in the anti-quark distribution with an 
emission of a boson. The term (III) is schematically de- 
picted a.sQ + b ^ q. This process has a support for w < 
and the quasi-quark Q is dominated by hole components 
of anti-quarks in this energy range. Through the process 
of the term (III), the quasi-quark is converted to the on- 
shell quark, and thus the mixing between a quark and a 
hole of anti-quark is induced. 

Scattering processes of a boson would induce a mixing 
between particle and hole states at finite T, which are 
called the resonant scattering [ll], [13, EBl • The resonant 
scattering causes a level repulsion in the energy spectra of 
the fermionic excitations. In the case of superconducting 
phase transition, for example, the preformed-pair modes 
induce a resonant scattering between particle and hole 
states in the Fermi sea; the outcome is a gap-like struc- 
ture around the Fermi surface, i.e., the so called pseudo- 
gap, in the fcrmion spectral function [1511. For the chiral 
phase transition, it was shown [Til . Il2| that a resonant 
scattering between a quark and an anti-quark hole is in- 
duced through the coupling with the chiral soft modes; 
there appears a three-peak structure in the quark spec- 
tral function p_|_ (w) in the low-energy and low- momentum 
region [ll|, see, Eqs. (3.18) and (3.19) in Ref. 



Now, let us consider the level mixing for the present 
case with m/ 7^ 0. We first focus on the term (II) with 
the energy conservation 

= ~Ef{~k) + Et{k) = c^(")(fc), (6.1) 

with k = |k| , where ujq is the energy of the quasi-quark Q 
with zero momentum, a'^^'^(fc) denotes the energy of the 
two-particle state composed of a boson with momentum 
k and an antiquark-hole with momentum — k, and is var- 
ied with k from M_ = rrib — m/ to as shown by the solid 
lines in Fig. [201 Although u;'^^^\k) can take any value be- 
tween and M_, there exists a value where the process 
corresponding to (II) is most probable. Such a value may 
be given by the peak position of Iml]^(cj(^^^ (fc)) due to 

the term (II). We denote the peak position as e'^^g and 
show it in Fig. 1201 Then, the physical energy spectrum at 
zero momentum is realized effectively through a level re- 
pulsion between s'^^g and the on-shell free-quark energy 

TO/. From Fig. [20l one finds that gets to exceed 
TOy at a mass between mf/mb = 0.2 and 0.3 when ruf 
is increased. This crossover is also seen from the move- 
ment of the peak position of ImI]^(Li;) relative to m^, 

as shown in Fig. [211 For mf /mi, = 0.1 and 0.2, e'^^g is 
larger than m/, while for mf/mb = 0.3, becomes 
smaller than mf. Now, when e'^^g > mf, which is the 
case for mf/mb < 0.21, the state with the unperturbed 
energy m f tends to have a lower energy than to f , while 

that with e^g-* to have higher energy than e^^-* by the 
coupling through the resonant scattering with the boson. 
On the other hand, when e'^^g < mf, as is the case for 
mf/mb > 0.21, the state with the unperturbed energy 
TO/ tends to have a higher energy than to/, while that 

with e^^g to have lower than e^'^g . Thus we have seen that 
the level-crossing phenomenon between poles (A) and (B) 
or the trajectory rearrangement at mf/mb — 0.21 can be 
nicely account for in terms of the level repulsion induced 
by the resonant scattering. 

Similar discussions hold for the process of the term 
(III) , where the energy conservation reads 

c^Q = S/(-k)-£;b(k)=c^(™)(A:). (6.2) 

^(iii)^^'j jg drawn in Fig. [201 with a solid line in the neg- 
ative energy region. The most probable energy, e^.™'* , in 

this case is then negative, and thus e'^g^^ < mf is satis- 
fied irrespective of the value of m/. The physical energy 
spectrum at zero momentum is again obtained effectively 
as a result of the level repulsion between s^^g^'^ and to/. 

Because e^™'' and to/ never cross, there is no qualitative 
change in the T-dependence of the pole (C) unlike the 
process of the term (II) : The real part of the position of 
the pole (C) monotonically decreases as r increases, as 
shown in Figs. [TTl [T5I and [T51 As mf/mb increases with 

T fixed, goes farther from to/, which means that 
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the effect of the level repulsion gets weaker and thus the 
strength of the pole (C) gets weaker, as shown in Fig.Ul 
We thus find that the T-dependence of the poles in 
the quark propagator for each ruf can be understood in 
terms of the level repulsion which is induced by a reso- 
nant scattering of the massive boson with a quark and 
an anti-quark hole. 



VII. SUMMARY AND DISCUSSION 

In this paper, we have investigated the spectral prop- 
erties of quarks with a Dirac mass to/ coupled with a 
scalar and pseudoscalar boson with a mass nib at finite 
temperature (T). We have employed Yukawa models to 
describe the system and studied the spectral function and 
the pole structure of the quark propagator for zero mo- 
mentum at one loop. The spectral function in the same 
models at finite T was studied for the case with to/ = in 
Ref. [12] , where the formation of the three-peak structure 
was observed at intermediate temperature T/rrib ~ 1. In 
the present study, we have examined effects of to/ on 
the spectral function in the range mf/rrib < 0.3 focus- 
ing on the spectral properties at zero momentum. It was 
found that, although the three-peak structure found in 
Ref. survives with finite to/ , it tends to be suppressed 
as uif /rrih becomes larger. We showed that the T depen- 
dence of the free pole at u; = to/ changes qualitatively as 
mf/rrib becomes larger. 

We have newly investigated the pole structure of the 
retarded quark propagator in the lower-half complex- 
energy plane and the residues. Our numerical calcula- 
tion has identified the three poles corresponding to the 
peaks in the spectral function, although near the imagi- 
nary axis, there are infinite number of poles which do not 
form peaks in the spectral function. The T-dependence 
of the three poles is consistent with that of the corre- 
sponding peaks in the spectral function. We have found 
that the T-dependence of the poles changes qualitatively 
around mf/nrn, ~ 0.21 (0.23) for Yukawa models with 
scalar (pseudoscalar) boson. This behavior is consistent 
with the behavior of the spectral function and suggests 
that the physical contents of the collective quark excita- 
tions at finite T change qualitatively around this value of 
the quark mass. We have given an interpretation for this 
characteristic change of the collective excitations in terms 
of the level crossing induced by the resonant scattering 

The spectral function in the Yukawa model at finite T 
was studied in Ref. [l^, for the case with TOf, = 0. In this 
case, the spectral function gradually changes between a 
single- and two-peak structures as T/TOf, is varied as an 
external parameter. It is an interesting project to explore 
the pole structure in this case, and the possible quahta- 
tive change of them in the whole region of parameters 
TO/, TOb and T. This result will be reported elsewhere. 

The present study has a relevance to the study of ef- 
fects of hadronic excitations [i, in the QGP phase 



on the quasi-particle picture of the quark. One of the 
plausible candidates of such a hadronic mode is the soft 
modes of the chiral transition Q since they can become 
light near Tc. The effect of the soft modes on the quark 
spectrum has been considered in Ref. 11] in the chiral ef- 
fective model in the chiral limit. In this case, the quarks 
are massless above Tc and the soft modes become mass- 
less at Tc. If the chiral symmetry is explicitly broken, 
on the other hand, the chiral phase transition at finite 
T becomes crossover and the quarks have a finite mass 
for any value of T [l^ . The soft modes do not become 
massless, either. These effects would suppress the forma- 
tion of the three-peak structure in the quark spectrum 
near Tc, as we have seen in the analysis of the present 
work. The masses of quarks and soft modes, however, 
varies dynamically as functions of T in the QGP phase. 
In order to investigate the effect of the soft modes more 
quantitatively, the study of the quark spectrum with a 
chiral model incorporating the explicit chiral symmetry 
breaking is desirable. 

The quark spectral function near but above Tc is re- 
cently analyzed in quenched lattice QCD in Ref. 
where the quark spectral function p+{oj) at zero momen- 
tum is analyzed as a function of ruf. It is shown that 
p+{uj) is well described by a two-pole ansatz. This result 
means that the three-peak structure is hardly formed in 
the quark spectrum in quenched QCD. We note, how- 
ever, that the effect of the hadronic excitations is not 
incorporated in the quenched approximation. The full 
lattice simulation would properly describe the effect of 
such excitations on the quark spectrum. 

It is an interesting problem to explore the effect of the 
multi-peak structure on experimental observables. For 
example, the dilepton production rate is an interesting 
candidate which can be affected by the structure of the 
quark propagator [iTj . The excitation spectrum at low 
energy may modify the thermodynamic quantities, such 
as the transport coefficients. To examine these observ- 
ables with the quark propagator obtained in the present 
study is beyond the scope of this work and left as a future 
project. 

Since the model which we have employed in this work 
is a simple but generic Yukawa model with a massive 
boson, the results can be applied to various systems at 
finite T composed of a massive fermion and a massive 
boson. For example, the quark spectrum might have 
a multi-peak structure through couplings with massive 
bosonic states instead of the soft modes, such as glue- 
balls, charmonia, and so on. Another example is the ex- 
citation spectra of neutrinos coupled with weak bosons 
at T near the masses of the weak bosons [3] • The multi- 
peak spectral function of a Dirac particle with a vanishing 
or small mass might be realized in some two-dimensional 
materials such as graphene 11911 and an organic conduc- 
tor a-(BEDT-TTF)2l3 salt which exhibit Dirac-type 
dispersion relations. 
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FIG. 21: The imaginary part of the quark self-energy ImE^(a;) for mf/mi, = 0.1,0.2 and 0.3 with fixed r — 0.5. The dotted 
lines represent lu — irif. r is chosen as the value at which p+{uj) and the position of pole (A) start to show a clear deviation 
from those for T = 0. 
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APPENDIX A: RENORMALIZATION OF THE 
T-INDEPENDENT PART OF THE QUARK 
PROPAGATOR 



In this appendix, we calculate the T-independent part 
of the quark self-energy I]j;^g(a;), which has an ultravio- 
let divergence. The renormalization is carried out using 
the subtracted dispersion relation. In the following, we 
concentrate on the quark part I]_|_.t=o('^), because the 
anti-quark part Ti^^T=oi^) be dealt with in the same 
way. 

We first note that ImE+^T=o is free from divergence 
and thus can be calculated without regularization. Tak- 
ing a limit T ^ O"*" in Eq. (|2.8p . one obtains 



subtracted dispersion relation 



ReI]+,T=o(^) =^ 



n 



-ci 



+ —P 



dC 



ImS 



R,T=0 



(C) 



(C-c.)(C-a)"' 
(A2) 



where a denotes the renormalization point, c; are the sub- 
traction constants, and n is the number of subtraction. 
To regularize the ultraviolet divergence, we use Eq. (|A2p 
with n = 2. 

The subtraction constants cq and ci may be deter- 
mined by the on-shell renormalization conditions for the 
mass and the wave function at a = m/, which are tanta- 
mount to requiring that 



and 



^+,T=oimf) = 0, 
(9S+,T=o('^) 



duj 



= 0, 



(A3) 



(A4) 



u!—7n f 



respectively. From Eqs. (|A3p and (IA4p . we obtain cq 
ci = 0, and then the renormalization is completed. 



APPENDIX B: GENERIC NATURE OF THE 
THREE-PEAK STRUCTURE IN THE QUARK 
SPECTRUM; AN ANALYTIC TOY MODEL 



Imi;+,T=o('^) 

(cj-Af_)(w + Af+) 



327r 



S{Lj)e^9{uj^ - Ml), (Al) 



where M_ = mi, — mf, = mi, + ruf, S{uj) = 

^ {oj^ — Af£)(a-'^ — M^), e^j = uj/ \uj\ is the sign of ui and 
0{uj) is the step function. 

To determine the real part ReS_(-.T=o('^)j we use the 



We have seen in the main text and Refs. [ll|, [12] that 
the two-peak structure in the imaginary part of the quark 
self-energy leads to the three-peak structure in the quark 
spectral function. In this appendix, we show that this has 
a generic nature using a simple but generic model. For 
simplicity, we consider the case for vanishing momentum 
as in the main text. 

First, let us consider the quark spectrum for the mass- 
less quark. To investigate effects of the peak structure in 
the imaginary part of the self-energy on the quark spec- 
trum, we assume that the imaginary part of the retarded 
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CO/oOn 



FIG. 22: The spectral function p+ and the self-energy E4 
7/ljo = 0.2. The quark mass is taken to be zero. 



for 



quark self-energy consists of two peaks expressed with 
Lorentzian forms, 



IniE+(w) = -M^ 



71 



72 



(Bl) 

where uJi^2 denote the peak positions and 71,2 denote the 
widths. An overall factor having mass dimension two 
is put because the quark self-energy has mass dimension 
one. Since we are interested in the number of peaks in the 
quark spectrum, the deviation of the peak shape from the 
Lorentzian form is not important here. For the massless 
quark, we have loq = lui = —lu2 and 7 = 71 = 72 from 
a symmetry property of the self-energy. As discussed in 
the main text and Refs. [l2|, the formation of the 
peaks in Iml]+ is caused by the Landau damping which 
is scattering processes with a boson. As T increases, the 
peaks become higher rapidly as shown in Figs. [5] and [6l 
which corresponds to 7/079 being smaller. 

The corresponding real part of the self-energy is 
uniquely determined from the imaginary part: Because 
the Landau damping is a medium effect, the real part is 
ultraviolet finite and given by the un-subtracted disper- 
sion relation. 



ReE+M^T'/ - 



doj' ImE_|_(a;') 



OJ' 



- Wo 



■UJQ 



_{UJ0 - LO)^ + 72 



{loq + ijjY + 7^ 



(B2) 



The spectral function is then obtained from the self- 
energy, pj^{ij}) = -I/tt • Iml]+(tj)/[(tj - ReE+(w))2 -f 
ImE+(w)2]. 

As an example, we plot the spectral function and the 
self-energy for 7/070 = 0.2 in Fig. [22l For simplicty, we 
have set M = ujq. From the number of crossing points 
between ReE_|_ and the line lo/ujo^ one finds that this self- 
energy is close to that for T = 1.5 shown in Fig. [51 We 



see that the resultant quark specral function p+ has a 
clear three-peak structure. For smaller values of 7/070, 
i.e. higher T, the three peaks are sharper and higher 
than those in Fig. [221 

For very high T, however, the strength of the central 
peak in gets weaker and eventually p+ becomes to 
have only two peaks which correspond to the normal 
quasi-quark and the anti-plasimino as shown in Fig.[31 In 
the Yukawa models, this can be understood from the fact 
that as T increases, in addition to two peaks in ImE-|_ be- 
ing higher, the positions of the two peaks approaches the 
origin. In terms of Eq. (|Bip . this means 07i/r, 072/T ^ 
for T 00, thus it becomes hard to distinguish the two 
peaks, and thus ImE_|_ only shows a one-peak structure 
as a whole. One can show in a generic way that the one- 
peak structure of ImE-|_ leads to a two-peak structure of 
the quark spectral function. In fact, this rule can be ex- 
tended to the cases where ImE-(- has multiple peaks; the 
number of the peaks in ImE+ is significantly reflected in 
the number of the peaks in p^. 




co/tOn 



FIG. 23: The spectral function p+ and the self-energy E4 
7/070 = 0.2 and mf/ioo = L5. 



for 



Next we consider the quark spectrum for the massive 
quark. In this case, 07i and 71 are in general not identical 
to —072 and 72, respectively. However, from Fig. [51 we 
see that such asymmetry is not so large for small values 
of the quark mass m/. Then, for simplicty, we employ 
the same self-energy as that for the massless quark. The 
mass effect only enters the expression of the spectral func- 
tion, p+{u!) = — l/vr • ImS+(o7)/[(o7 — my — ReE4.(o7))^ -t- 
ImE+(w)2]. 

As an example, we plot the spectral function and the 
self-energy for 7/070 = 0.2 and nif/ujo = 1.5 in Fig. [23l 
From the number of crossing points between RcE+ and 
the hne {lu — mf)/uio, this self-energy is colse to that 
for T = 1.5 in Fig. [6l We see that the rightmost peak 
is enhanced, which results from the shift of the crossing 
points. This result of the quark spectrum is also consis- 
tent with that in Fig. [H 
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For very high T, the effect of the quark mass be- 
comes relatively smaller and thus the quark spectrum 
approaches that for the massless quark mentioned above. 

In short, we have shown that a simple but plausible as- 



sumption on the peak structure in the imaginary part of 
the self-energy generically leads to the multi-peak struc- 
ture in the quark spectral function obtained in dynamical 
models such as the Yukawa models. 
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